1 



m 
o 
o 

(N 

Oh' 
(D 

in 

(N 



> 

oo 
(N 

a^ 
o 

m 

o 



X 



Transverse Momentum Distribution in 6 ^ 57 
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We present the complete calculation of the transverse momentum distribution for the decay b 37. The 
contributions of the leading operator O7 are computed: infrared logarithms are resummed with next-to-leading 
accuracy according to usual techniques of resummation. Non logarithmic terms are evaluated to 0{as) by 
calculating one loop diagrams. 



1. Introduction 

The rare decay b —^ sj has been widely studied 
in the last years, both experimentally and theo- 
retically, for several reasons: being loop mediated 
even at the lowest order, it was used to search 
signals of new physics. Moreover it is a way to 
measure CKM matrix elements, which is one of 
the main goals of beauty physics to constraint the 
Standard Model. 

Pcrturbative QCD has been applied to b physics 
since a long time, because the QCD coupling con- 
stant is assumed to be small at b quark mass 



as(m2)^ 0.21 



(1) 



However a direct check of pcrturbative QCD is 
difficult to perform, because several poorly known 
parameters are involved in theoretical predic- 
tions (for example CKM matrix elements and the 
beauty quark mass) . On the contrary one usually 
prefers to assume that pcrturbative QCD is reli- 
able in b physics and can be used to extract these 
parameters. 

In our opinion |T] it is not clear whether the pcr- 
turbative approach is viable at 6 mass energies, 
both in inclusive quantities and, a fortiori, in 
semi-inclusive distributions. 

In the latter case a powerful tool is the resum- 
mation of the large logarithms appearing at the 
border of phase space. This approach had a wide 
application in past years in processes such as 
e+e^, pp, deep inelastic scattering, . . . , for semi- 



inclusive distributions and in recent years has be- 
come an ordinary technique in b physics too. 
It is well known that the resummation of large 
logarithms for a generic distribution D{x) gives 
the general result [2]: 



D{x) = K{as)^{x; as) + R{x) 
where 



(2) 



• T,(x; as) is a universal, process independent 
function resumming the infrared logarithms 
in exponentiated form. In particular 

logS(a;;as) = Lgi{asL)+g2{asL) + . . . (3) 

being L a large logarithm. The functions gi 
have an expansion as 



(4) 



and resum logarithms of the same size: in 
particular gi resums leading logs log"^^ 
and 52 the next-to-leading ones ag log" ; 

• K{as) is the coefficient function, a process 
dependent function which can be calculated 
in perturbation theory by evaluating the 
constant terms at the required order 



K{as) = 1 



Cpas 



fci + 



(5) 
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• R{x), the remainder function, is process de- 
pendent and satisfies the condition 



R{x) ^0 for x^O 



(6) 



It takes into account hard contributions and 
can be calculated in perturbation theory as 
an as expansion. 



R{x) 



asCf 



(7) 



For the previous reasons, a precise calculation of 
the terms involved in this approach is needed to 
compare theoretical predictions and experimental 
data. 

For b s"f the resummation has been applied 
in the photon spectrum distribution which is 
sensitive to longitudinal degrees of freedom. We 
have extended the analysis to the transverse mo- 
mentum distribution of the strange quark with 
respect to the photon direction, which is sensi- 
tive to transverse degrees of freedom 4 . 
The resummation of large logarithms seems to be 
required in the transverse momentum distribu- 
tion for b ^ sj, because the double logarithmic 
correction becomes large near the border of phase 
space, compared to 1 (lowest order contribution 
to the rate): 



1 Cpasiml) 



log' 



A2 



-0.7 



(8) 



if we push the transverse momentum to its phy- 
sical limit, A ^ SOOMeV, before the appearance 
of Landau pole effects. 

The single logarithm becomes large too, due to 
its large numerical coefficient 



4 n 



0.6 



(9) 



However hard contributions as the function R(x) 
should be carefully evaluated, being at these ener- 
gies more relevant than in processes such as e+e~ , 
deep inelastic scattering, due to the size of ^]. 

2. Transverse Momentum Distribution for 

b ^ sj 

Let us consider the decay b ~f sj in the b 
rest frame and define the photon direction as the 



z — axis. Let us consider the transverse momen- 
tum p± of the strange quark with respect to the 
photon direction and, in particular, the partially 
integrated distribution 



r 1 , , 



(10) 



where x is an adimensional variable defined as 

(11) 



Pi 



being rri},, the b quark mass, the hard scale of the 
process. Fq is the Born amplitude, as defined in 
0. 

According to an accurate calculation of D{x) 
requires two steps: 

• the resummation of large logarithms aris- 
ing at the border of phase space (infrared 
region) 

• a fixed order calculation to extract non lo- 
garithmic contributions, such as constants, 
contained in K(as), and functions vanish- 
ing for a; ^ 0, contained in R{x) (J. 

2.1. Resummation of infrared logarithms 

The resummation of infrared logarithms with 
next-to-leading accuracy requires the knowledge 
of the coefficients of double and single logarithms 
at one loop and the coefficient of leading loga- 
rithms at two loops. 

At one loop the distribution takes the look 



D{x) = 1 — as^ilog^a; -I- as -Bi log a; 
+ ^k, + ^r,{x) 



(12) 



The coefficients Ai and Bi can be calculated by 
general properties of QCD, applying the pertur- 
bative evolution for the light quark, described by 
the Altarelli-Parisi kernel, and the eikonal ap- 
proximation for the massive quark 
The coefficients turns out to be U 



TT 



4 TT 



(13) 



Let us notice that the coefficient Bi is larger than 
in processes involving light quarks only, due to 
gluon emissions from the heavy quark and this 



3 



produces an enhancement of the single logarith- 
mic term. The two loop coefficient A2 is known 
from literature \7\ 



C'f, ,67 



9 



(14) 



The resummation of transverse momentum is per- 
formed in the space of impact parameter b 



dp^e^v-.-i (p^) (15) 

1 dp±_ 



which allow the factorization of the phase space 
and of the kinematical constraint. 
The result for the functions gi and 52 appearing 
in (PJ is 



52 (w) 



^-[log(l 
2/3o oj 



(16) 



2PI 



1-UJ 

log(l 

l-LO 



+ log(l-w) 



1-OJ 



+ llog2(l-c.)] + ^log(l-c.) (17) 

being u = /3oaslog(^) and 60 = ie''^''"' . 
"0 

Let us notice that the resummed distribution di- 
verges for w — > 1 and this signals the appearance 
of non perturbative effects, related to the Landau 
pole of the theory. They appear when the distri- 
bution approaches the limit of the phase space in 
the infrared region, that is, roughly speaking, for 
Pjl A. 

The resummed distribution is more reliable than 
the fixed order calculation near the border of 
phase space, nevertheless it breaks down for small 
value of the transverse momentum. 
This divergence should be factorized in a non per- 
turbative structure function: for the photon spec- 
trum an approach based on an effective theory 
has been used, by introducing a shape function 
[HI, but in this case we conclude in 0] that this 
approach is not viable and one should try another 
way to include non perturbative effects. 

2.2. Fixed order calculation 

A complete resummed distribution with next- 
to-leading accuracy requires the calculation of the 



constant ki in (|12|1 . The first term in the expan- 
sion of the remainder function, ri{x), should be 
relevant too, due to the low value of the hard 
scale. They can be extracted by an explicit cal- 
culation of real and virtual diagrams at one loop. 
The process 6 — > S7 is loop mediated in the Stan- 
dard Model: an effective hamiltonian has been 
developed to integrate out the heavy particles in 
the loop |S] 



G ^ 



(18) 



where Cj are Wilson coefficients and Oj{ii) are 
effective operators (i = 1, . . . , 8) 9 . 
O-! is the most relevant operator because it is 
affected by a logarithmic enhancement in the 
infrared region, while the other operators con- 
tribute to constants [HI] . 

Real diagrams are calculated in dimensional re- 
gularization in n — 4 -I- e dimensions. 

Db{x)^ I d<^ MA+e{uJ,t) 6[x-uj'^t{l-t)]{l%) 
Jo 

being $ the 3-body phase space, u the gluon en- 
ergy fraction and t = {1 — cos0)/2, where 9 is 
the angle between the gluon and the direction 
— z. A^4+c(tJ,t) is the matrix element with the 
insertion of Oj in the vertex, in dimensional re- 
gularization. 

The integration may be analytically performed by 
introducing harmonic poly logarithms as in |11| . 
satisfying the properties 



V 



dy' g{a; y') ay^O 





logy 

V 



dy' 9{a;y') H{m^^i;y') 







—H{m^;y) = g{a;y)H{m^^i;y) (20) 
dy 

where the basis of functions g{a; y) is 

y 



y+l 



4 



5[-3;?;] 



2(1- 



(21) 



Details of the calculation and the explicit result 
for the first term in the expansion of the remain- 
der function are shown in 

Virtual diagrams are calculated with ordinary 
techniques: they cancel the infrared poles ari- 
sing from real diagrams and don't depend on the 
kinematics, giving contribution to the coefficient 
function only. In particular the vertex correction 
is calculated by reducing the diagram to simpler 
topologies and using integration by parts identi- 
ties: details of the calculation are discussed in pp. 
Summing real and virtual diagrams the first term 
of the perturbative expansion of the coefficient 
function for the operator Oy turns out to be 



log ) 



(22) 



An improved result for the coefficient function, 
taking into account the other operators of the ba- 
sis (|18|) . can be found in 

3. Conclusions and Outlook 

h — > S7, though a rare decay, seems to be a 
very clean process to study hadronic physics in 
h decays and to apply resummation techniques 
to check whether they are reliable in h physics 
or not. A very accurate calculation of all the 
ingredients involved in the resummation " recipe" 
is needed to this aim. Otherwise, it would be im- 
possible to conclude whether discrepancies from 
experimental data may be attributed to this ap- 
proach or to other sources of uncertainty typical 
of b physics, such as an imperfect knowledge of 
CKM matrix elements or the h quark mass. 
A more interesting application should concern 
the transition b c, which is phenomenolo- 
gically more relevant, but whose calculation is 
complicated by the presence of the c quark mass. 
From this point of view the study of 6 —s- 57 
may be seen as a first step to apply resummation 
techniques to b physics, before going to more 
complicated but more important processes. 



The contents of this paper arise from my work 
with Prof. L.Trentadue and Dr. U.Aglietti: I 
wish to thank Prof. Trentadue for his collabora- 
tion and Dr. Aglietti for his essential help. 
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